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THE ANALYST. 

Vol. II. Nov., 1875. No. 6. 

ON THE DEVELOPMENT OF THE PERTUBBATIVE 
FUNCTION IN PERIODIC SERIES. 



BY a. W. HILL, NYACK TUENPIKE, N. Y. 

1. There are two modes of developing this function. In one, the 
numerical values of the elements involved are employed from the outset, and 
the results obtained belong only to the special case treated. This mode has 
been, almost exclusively, followed by Hansen, and is, perhaps, to be recom- 
mended when numerical results are chiefly desired. In the other, all the 
elements are left indeterminate, and thus is obtained a literal development 
possessing as much generality as possible. Certain investigations, arising 
from Jacobi's treatment of dynamical equations and Delaunay's method in 
the lunar theory, have invested the latter mode of development with addi- 
tional interest, and with it we shall be exclusively engaged in this article. 

In Liouviile's Journal for 1860, M. Puiseux has given us two memoirs 
on this subject, in which appears the general term of this function, but his 
formulae seem susceptible of modifications which would render them much 
simpler. More recently, in the volume of the same journal for 1873, M. 
Bourget has presented the development in a more concise form by employ- 
ing the Besselian functions, but as he discards the use of the functions 6, (i) , his 
formulae, on this account are more complex. It is hoped, that, even if the 
expressions, given hereafter, are deemed too cumbrous for practical use, they 
may still possess some interest from a theoretical point of view. 

2. It is known that if we have a function S of a variable C, which is 
never infinite, and such that the relation 

function (£ + 2wr) ~ function (£) 
is satisfied for all integral values of i both positive and negative, it can be 
developed in a series of the form 

I t . (Al c) cos iC + Iv< s) sin iQ, 
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in which i denotes a positive integer; and that, in the eases where this series 
is infinite, it is convergent. 

In general the handling of periodic series is easier if we introduce imag- 
inary exponentials in the place of the circular functions. Thus, e denoting 

the base of natural logarithms, we shall put 2 = e , whence 

2 cos £ = 2 -f- z" J , 2 cos if — z l -+■ z~ l , 

2 1 /(— 1) sin C = 2 — z" 1 , 2 T /(— 1) sin i£ = z* — 2-* , 

2 = cosf -f l/( — 1 ) sin <£, 2* = cos if -f- j/( — 1) sin if. 

The above theorem then comes to the same thing as to say that S is de- 
velopable in a series of the form 

where the summation is extended to negative as well as positive values of i. 
The coefficients K are given in terms of the coefficients G by the equations 

kp = q + c_,, 

KP = (Q - C-V-l, 
except the case where i = 0, when K§ = C . It will be seen that when 
<S is real, Q is a complex number a + 6]/ — 1, and CLj, its conjugate a — 
b]/ — 1, which renders the coefficients JTreal, as they should be. 
The integral 

JVdf = /( cos if + t/(— 1) sin »c)dC» 

taken between the limits and 27r, vanishes in all cases except when i — 0, 
when its value is 2tt. Hence any function, capable of expansion in a series 
of positive and negative integral powers of 2, integrated with respect to f 
between these limits, gives, as the result, 2n times the coefficient oi :'s° in its 
expansion. And as the coefficient of z° in the function Sz~' is evidently 
Q, we have 



1 /*2ir 



This equation holds for all values of i, negative as well as positive, zero 
included. 

3. Let us now suppose that f denotes the mean anomaly of a planet, 
and let u be the eccentric anomaly, connected with the former by the equa- 
tion, e being the eccentricity, 

u — e sin u = f . 
In like manner as for f, we will introduce the imaginary exponential 

s = e , thus, as the last equation can be written 
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(u — eamu)i/ — 1 ty— 1 

£ = S , 



by the introduction of the variables s and z, this becomes 



SS 2V »' = 2, 



which is the transcendental equation connecting s and z. We have 
dC, = (1 — e cos «)dtt r= 1 — |(s + -1 \du. 

Substituting these values in the equation giving the value of Q, and no- 
ticing that, as £ and u both take the values and 2tc together, the limits of 
integration, when u is the independent variable, are the same as for £, we get 

But, from what precedes, we conclude that the coefficient of s* in the expan- 
sion of any function W, according to positive and negative powers of s, is 

1 r in 

±- { Ws-'du. 
2nJ 

Thus, from the foregoing expression for Q, we derive the following proposi- 
tion : — 

i being a positive or negative integer or zero, the coefficient of z', in the de- 
velopment of S, according to the powers of z, is equal to that of s 1 in the de- 
velopment of 



s ^)Q_. ( . + lj], 



according to the powers of s. 

As most of the functions 8, which are presented by astronomy for develop- 
ment in powers of z, are quite readily expanded in powers of s, this theorem 
is of much use. Another form can be given to it. For we have, integra- 
ting by parts 



jSz-'dc = — j/(— 1) J#z- ( ' +1 



\k 



= V==±Sz-<-^±C d fz-*dz. 
i i J az 

Taking the integrals between the limits £ = and £ = 2jt, we get 

t/-i r dS ~-^ 



r — — y— l C 

2ircJ 



- z 



2inJ o 
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Whence we conclude this proposition : — 

The coefficient of z 1 in the development of S according to the powers of z is 
equal to that of ' s 1-1 in the development of 

\MJiis-r) 

i ' ds 
according to the. powers of s. 

This theorem however is not applicable when i = 0. 

4. We shall often have occasion for the expansion of the function 

in powers of s; let us, for simplicity, put / = \f, and 
A(s-A) x s _^ v Tm ; 

e v « ' = S .S ■;=r2 i .J W 8'. 



We have 

e' As . e- 



, ... ^i ^ _ + __ + __ + . . . j 

L s 1.2 s 2 1.2.3 s 3 " r ' ' " J' 



X 
whence we conclude that 

jc>= _^_ n - -*- + *- - 1 

X 1.2...; L l-(*+l) 1.2(t+l)(t+2) "•J" 

This series is not applicable when i is negative; but if, in the function 

> k 1 

s "* • £~~, we substitute - for s, and change the sign of /, the function re- 
mains unchanged, hence 



and consequently 



X -I 



X -x x ' X ' 



by which the values of these functions for negative values of i can be de- 
rived from those in which i is positive. These functions are known as the 
.Besselian. By putting 

2J = 1 » • * 



l.(t+l) 1.2(«+l)(*+2) ' ' 
one will have no difficulty in deducing the equation 

5. We come now to the more complex function S of two variables J 
and C'; it is known that when this is never infinite and is such that 
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function (£ + 2wr, £' + 2i'?r) == function (£, C') 
it can be developed in a series of the form 

2V „ [Kg, cos « + t'O + -S^-i- sin « + i':')], 
where to one of the quantities i and i', we need assign only positive integral 
values, but to the other both positive and negative values. If we adopt 

another imaginary exponential z' — e , this is the same as saying that 

Cf — V fi 7 i,li< 

where the summation is extended to all integral values positive and nega- 
tive for i and V . Since we have 

2 V' ; ' = (cosiC + y (— 1) sin iC) (cos %%' + \/{— l)smi'?) 
= cos (ir + i'r>) + v /(— 1) sin {iC + i'f), 
the relations, which connect the coefficients if with the coefficients C, are 
-tvi'ji = Cii, -f- C _,-_(/ , 
*& = «?,,, - (?_,,_,, )t/— l, 

unless i and -i' are both zero, when 

JL o, o — °o, o • 
A course of reasoning, similar to that in the case of one variable, establishes 
that 

which holds for all integral values of i and i', positive, negative and zero. 

6. Supposing that £' denotes the mean anomaly of a second planet, 
whose eccentricity and eccentric anomaly are respectively e' and u', we have 

u' — e' sin u' = C.' , 

and by the adoption of the imaginary exponential s' — e , this is 

transformed into 

It is not difficult to see that we have the following theorem : — 

The coefficient of z'z"' in the development of S, according to the powers of 
z and. z', is equal to that o/s's"' m £/i,e development of 

according to the powers of s and s'. 

7. After these preliminaries relative to the general development of 
functions in periodic series, we come to the matter more immediately enga- 
ging our attention. The perturbative function for the action of a planet, 
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whose mass is m', on another, whose mass is m, is usually written 

, fl _ r cos f ~\ 



R 



R x — m 



'1 r' cos ^~| 
A r 2 J' 



where // denotes their mutual distance, ^ their angular distance as seen from 
the sun, and r and r' their radii vectores. The problem proposed is then 
to develop these two functions in series whose general term is of the form 
C it viz' 1 '. To this end it seems better to discuss the two portions of the 

general perturbative function, — and — Tz — , separately, and not, as 

most investigators, attempt, by a particular notation, to combine, in a whole, 

these two parts. Thus, in developing -, we shall have the term common 

to both functions, and may suppose that r' denotes the radius vector which 
belongs to the planet more distant from the sun. But, in treating the sec- 
ond part, we shall suppose that r' belongs to the disturbing planet. The 
following equations are well known, 
j2 __ r /2 — 2rr' cos <p + r 2 , 

cos ip — cos (v + 77) cos (v' -f- H) -f- cos 1 sin (v + II) sin (v' -\- II) 
= cos (v — v' + 77 — 77') — 2 sin 2 1 Jsin (v + 77) sin (v' + II'), 
where v and v' are the true anomalies, and 77 and 77' are the angular dis- 
tances of the perihelia from either point of intersection of the planes of the 
orbits, and I is their mutual inclination. 

8. Attending then, in the first place, to the development of -,, we have 

to notice what are the conditions under which this quantity can be developed 
in powers of z and z'. In the case of two elliptic orbits, the only one we 

shall consider here, it is plain that -, is always finite and continuous provi- 
ded the orbits have no point in common. Here we must make two cases 
according as the value of sin I is not or is zero. In the first case it is evi- 
dent that the orbits can meet only on the line of intersection of their planes. 
Hence, p and p' denoting their ssmi-parameters, there will be two, one or 
no points in common, according as two, one or none of the equations, 
p'(l + e' cos 77')- 1 = p(l + e cos 77)- 1 , 
p'(l — e' cos 77') -1 = p (1 — e cos 77) _1 , 
are satisfied. In the second case, where the orbits lie in the same plane, 
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there will be two intersections or none, according as the equation 
p'[l + e'cos(/ — w')] _1 — j°[l + ecos(A — «*)] -1 , 
k being the unknown quantity and to and to' the longitudes of the perihelia, 
admits real or imaginary roots. If we put 

pe' cos to' — p'e cos to = A cos a, 
pe' sin to' — p'e sin to = A sin a, 
this equation takes the form 

A cos (A — a) = p' — p. 
The roots of this are imaginary when 

(p' — pf > p 2 e' 2 — Ipp'ee' cos (to — to') -f- jj'V. 
9. If we put 

P = r' 2 — 2rr' cos (w — »' + 77 — 77') + r 2 , 
§ = 4 sin 2 J J . r sin (» + 77) . r> sin (»' + 77'), 



we have 



^ 2 = Prft 

4 = cp + « _i 



j 



= P"2 — 1 P~§Q 4- I § P~^0 2 
2 "C 2"4 K 



a series we shall denote thus 

-.yj 1 ■ O ■ . . ^K X ) p 

I O A 07. 



1 _ *5W_iv 1.3...(2ft-l) p -^±l 



z/ A, = 2.4... 2A 

10. In order that this development of - in a series of ascending pow- 

ers of §, or, if one likes, of sin 2 \I, may be legitimate, it is necessary that 

the elements of the orbits should be such that the numerical value of -~ 

should be always less than unity. P is the square of the distance of the two 
planets after the plane of the orbit of one has been brought into coincidence 
with the plane of the other by revolving it about the line of intersection of 
the two planes. Taking then a system of rectangular axes passing through 
the centre of the sun, and directing the axis of x along the line of intersec- 
tion, it is plain the equations of the orbits may be written 

y (x 2 + f) + ax + @y ~ p, 

V(x' 2 +y' 2 )+a'x' + P'y'=p', 

a, /S, «', )9' being constants. And the variables x, y, x', y' satisfying these 
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equations, the question depends on the finding of the values of them which 
render the expression 

d = yyL 

yx — x'Y -t- (y — y'f 
a maximum or a minimum. According to the known theory of maxima 
and minima, the equations, which, in combination with the equations of the 
orbits, give these values, are 

x . 



= o, 



y>- 2D(y-y')+p [_^i ? . + ^] =,0, 

where ji and // are the multipliers of the partial derivatives of the two equa- 
tions of condition. A complete investigation of this question would be 
conducted in the following manner. Eliminate from the seven equations 
last given the six quantities x, y, x', y' , p, // ; the result will be an algebra- 
ical equation determining the unknown D Having derived the Sturmian 
functions of this, one will ascertain by the substitution of the values _D = 

D = + oo, and again of D — r-^~JTr> D= — o>o, whether 



4 sin 3 !/' ' & 4sin 3 JI 

any roots lie between these limits ; if none, — can be expanded in a series of 

ascending powers of sin 3 \I, in the contrary case not. In this way we shall 
arrive at the condition or conditions necessary and sufficient for the legiti- 
macy of this expansion. 

11. This procedure would doubtless lead to very complicated formu- 
lae, hence we are obliged to pass over it. However equations can be read- 
ily got, which, by a tentative process, afford the maxima and minima values 
of D. Multiply the four equations last given respectively by x, x', y, y' 
and add the resulting equations, having regard to the equations of the orbits 
and the value of I), we thus arrive at the simple relation 

pp -+- p'p' — 0. 
Putting, for simplicity, x — r cos 6, x' = r' cos 6', the addition of the first 
and second of the same group of four equations gives 

p (cos d + a) + p' (cos 0' + a') = 0. 
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By combining this with the preceding is obtained 

cos tt' + u! cos d ~f « 

P~' P 

Again the addition of the same equations, multiplied severally by x, 
— x', y, — y', gives the equation 

2D (r' 2 — r 2 ) = p'ft' — p/u. 
Dividing the left member of this by 2D (x' — x), and the terms of the right 
member by its equivalents derived from the first and second equations, we get 

r' 2 — r 2 p' p 

T> i 7,i t~ 7.7Z7J \~~~7J 



x' — x cos 8' + a' cos -j- a ' 



r'cosO' — rcosO __ cos -f- a rp,. 



and the equation 



r' 2 — r 2 2p 

cos If -f- a' cos -f- a 



P P 

determine the values of the variables 6 and 0' which render D a maximum 

or minimum. When the orbits are nearly circular these values are in the 

neighborhood of \n or \x. When both orbits are circles the solution is 

very simple, and we find that in order the development may be legitimate, 

we must have 

. I . n! — a. 
sin - < 



2 " 2 l /{aa')' 
a and a' being the mean distances of the planets from the sun. 

12. Assuming that this development is legitimate, we have to develop 

P 2 Q h in terms of s and s'. We have 

r cos v = a (cos u — e) — -^(s + - — 2e), 

Z\ S I 

r sin v = rtj/(l — e 2 )sinw = j/(l — e z )U — -j, 

whence 

r cos •» + r sin « y 7 — 1 = rs 

and by putting 

2 "*' l + T/(l-e 3 ) <"' 

we get re = ays 1 1 — - ) . 
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And the value of ?•£ is evidently obtained by substituting in this - 

8 

for s, hence rs l = arj - (1 — om?. 

s 

From these two equations may be derived 

vy'—l __ g — a) 
1 — cos 
Writing y for // — //', we have 

...2*+l 

( r /-2Pj -5— _ 1 j __ 2 il C os(»— »' 4- r) + -Ti 

Tlie right member of this is developable in a series of integral powers of the 
exponentials ' when -£- is always less than unity. This con- 

dition is fulfilled when we have a(l + e) < a' (I — e'). Writing g for 



e , let 



(./) 



_2&+l i=+<*> <i> -ii'-i.— i/+y)i/— 1 

J"=-c» 2 






U> 



Bp s+ i is the same function of ^— that Laplace's 6 2 s+1 isof— 7 = a. The 
approximate value of — being a, any function of — T can be expanded in a 



series of ascending powers of ~ — a by Taylor's Theorem. And as we have 

r(l — ws)(l — — ) 1 

V Ll — ll, 



r 

77 



a — a 



consequently 



d n b 



(j) 



— j*- »= n! da" 



2&+1 
2— -{ 



i? (1 — <«8 )(l 



,'(l_a.'«0(l— y 



1 
1 \, 

I 
1 

J 



n being an integer, and n\ denoting the product of all integers up to n in- 
clusive, it being understood that 0! = 1. Expanding the last factor of this 
expression by the binomial theorem, and employing the notation [i, j] for 
the coefficient of x s in the expansion of (1 -f xf, we have, p being an integer, 
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13. In the next place the development of Q in terms of s and s' must 
be formed. We have 

. . m 1 f (« -h n )y--i -(« + n >,/-] 



2j/ — 1 



r'e —re , 



and putting // + //' = 6, h = £ , 

we find 

§ == _ aaV/' sin 3 ^. |j» (l - -^j — — (l — wsj V 1 ^ 1 ! 

^E^-t)'*--?! 1 — ¥ ) s " 

Raising this expression to the ft th power, and multiplying by 

r (2*+i) r- , «A-]-<2*+i) 

r' = ay(l — a>V)Il — — J , 



we 



find that the part of r'- (2 * +1) §* which has A*'" as a factor is 



y.y/Y-w-i) s i n 2^ v (— lf"[ft, «'][/«, * — *'" — «'] 



Ol 



\ 2»»'+2n' 



1 T jfcS-iJII 

A,aW-< ; ' +1 >sin 2 *4 - 
a' ' 2 .,.-=o 

X s 2'"'-*+2»' (i __ B,«)a»-s«"'-a«'/i _.. 

14. We are now in the possession of all the developments necessary 
for exhibiting the function - in terms of s and s'. In order to obtain the 

part of this function which has g l "h 1 '" for a factor, we must put, in the for- 
mula of § 12, 

j ~ i" __ i'ii + ft _ 2n', 

and the chief operation here is the addition of the exponents of the quan- 
tities s, 1 — cos, 1 — — , and the similar functions of s' which are found 

s 

_2*+i_ O) 
in the three formulae for (r'-~ 2 P) 2 , B 2h+1 and r'~ (2 * +1 > Q k . For 

2 

brevity we will write 

Li,J 2.4 2ft " 
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Then the part of - , which has g^'h 1 "' for a factor, is 

1 fa* n,=k-im n=* p=„ k-iW+n-p g.^ ,,/j ^ fc _ j//> — n']j>, j)] 

2a' *•=>"» n'=o «=o p=o n! 

X** w . sin-- 

w *+J-H"+i'" 



"A" 






We observe that in this expression the summation with respect to n' af- 
fects only the integral coefficients [k, n'~\, [k, k — %'" — n'~] and the upper 
index of the quantity b, hence if a new function of a is assumed, which is a 
linear function of the b's, and such that 

(i",i>") nl=k-itll (i"~i'"+*- 2»0 

Bm = 2 [*,<] [M-* v " -n']6 a n-i 

— 2— »'=0 — 2~ 

it will take the following simpler form 

— 2 2 2 ( 1) 1>J L" > Pi 

2a' k-=wt n=o P =o n! 

(»", i'") 

X«'' +n -rir- 2 sin 2 * i 
aa" 2 

Xy /+^ s "'+""(l — ««)*+*-<"-'" M — —\ 

X 7 '-*-*-i 8 '-"-H'"(i_ w V)-*- J, - 1+i "-""(l — ^) ^"tf"'. 

15. In order to get the coefficient of s' z /,f in the expansion of -, ac- 
cording to the foregoing investigation we must multiply the preceding ex- 
pression by — 7 e 2 v s ' ' 2 v «". 

Hence, if for brevity, we adopt the functional notation 



s Q ^^(l-o^y^i-v) 

£' f j ) = 7' V (1 -o-V^fl - ^ 



.2 



A «"+J *« 



<-|) 



("4) 



the coefficient of z* z* ^ { " h 1 '" in -, will be equal to the coefficient of s i s' v in 
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1_ 2'°° V" 2'" ( i\!'-i'"+t~p lALyVRJ 

2a' iZziin 7=0 p=o nl 



d n R (i "' "'"^ 

Xa*" 1 "" L 2 * +1 sin 

d,a n 2 



/ / k+ P \ /-{k+p+l)\ 
2i ~ . Si i"+i"' J .S'( -{i"-V») J . 



If then the coefficient of s* in the expansion of 8 is denoted by £ follow- 
ed by the same indices, and the coefficient of s"' in the expansion of <S" by 
£' in like manner, £ will be a function of e only, and £' a function of e! 
only; and, it being understood that each argument is taken but once, that 
is the negative of the argument is not considered, the coefficient of 

cos ( i£ + i'C' •+ i"r + i'"0 ) 
in the expansion of -. is expressed thus 



(-If 

d nDd»,i»0 J I ic+p \ /-(i+p+1) 
** n z-_l_ 1 • o I- r" I I n I 



a' kzznii »=o pzzo nl 



X«* + "-^-?i±± sin 2 * 2 



/ k+p \ ,-(k+ P +l)\ 



da" 2 

As in this formula h ought to be a positive integer, it will prevent em- 
barrassment, if the arguments are so taken that V" may not be negative. 
In the case where i, i', i" and V" are all zero, the expression must be divi- 
ded by 2. 

16. Thus we have arrived at an expression for the general coefficient 
involving only three signs of summation; and it may be remarked that all 
the coefficients are exhibited in precisely similar forms. Thus, to pass from 
one argument to an other, we have only to make the suitable changes in the 
two lower indices of the functions £ and £' and in the upper indices of ff, 
and commence the summation with reference to k with the new value ofi'" 
instead of the old. Hence, from this expression, we can write out a scheme 
or blank form, which, when the indices proper to the argument are filled 

in, will be the coefficient of the cosine of it in the expansion of -.. Such a 

blank form is written below; the indices i" and %'" are omitted from /?, 
and the two lower indices from £ and £', and the upper indices of these 
quantities, for the sake of facility in writing, are placed to the right and at 

the foot. The factor — , common to the whole expresson, is also omitted, so 

a' 
that the formula gives the coefficient in the expansion of -p In making use 

of it, one must commence at the portion which has sin 2i '" J7 for a factor, all 
the preceding parts being supposed to be suppressed. As to the double sign 
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before these factors, the upper or lower is taken according as i'" is even or 
odd. It is hoped that a sufficient number of terms have been written to 
render the law evident, so that they may be continued as far as desired. 

K£ F ' 



1 da 



._ o -1 



■£ F 

1 -2 . 



+h>°^[EF 2 E F + E F 

1.4 aa {_ o -l 1-2 2 - 



1 sd^j 
"1.2.3 ° da* 5 



££' -3£F +3££' -££' 

L. -1 1-2 2-3 S -4 



1 • , J 

2 Sm 2 



2 1-2 



1 ft2 d|? 
1 da 



£ F -EF 

.1.-2 2 -3._J 



-r- 9 « 3 %r^^ - 2 ££' + ££"' 

l.J da |_ 1-2 2-3 3 - 



T^4!*[ £ i £ - 2 -- 3f 



£ +3££ -££ | 



3 -4 



+ 



.1.3 • t I 

~i~ sin — ■ 

~2.4 2 



c?BsE F 

? 2 -3 






^1.2 da 2 



1 J>B 
1.2.3 da 3 



v 1 



'iYe? -3£F +3££ -££ J 

l_ 2 —3 3—4 4—6 5 — 1> _J 



+ 



1.3.5 • 6 7 
: 2X6 Sm 2 



0?BnE F 

2 3-4 

\^\\£ F -E F 1 

1 aa 3-4 4-5—1 



+ a « 5 %r^ ^ - ^ ^ + e f i 

1.2 aa'' L 3 -4 4-5 5 -e_l 

-l**M E p _ 3 f£' +3££ -££' I 
1.2.3 da? 2 L 8-* 4 -« 5 " 6 s ~ 7 ~ J 
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For illustration, let it be desired to obtain the coefficient of 

cos (2£ — b£> -f 2f), 
from which arises the larger part of the great inequality of Jupiter and Sa- 
turn, we have onlv to imagine that the lower indices [ 2 are every where 



B) 



ndices [ 2 J 



applied to £, and the indices f — 2 J to £' , the indices (2, 0) to B> 

and as we have V" — 0, we suppress nothing and take the upper of the 
double signs. 

17. The quantities B are very simply expressed in terms of the 6's, 
the following are all that are needed when terms of the eighth order with 
respect to the inclination of the orbits are neglected. 

i, 0) (>') 

Bx = b v 

(•', 0) (>+l) (i-l) 

Bz = h + b z , 
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2 2 

18. In computing the factors of the preceding formula which depend 
on £ and £', the following abbreviation can be used. M„ denoting the fac- 
tor which multiplies = -\ a k + n d "fy k +\ 

n\ da L 

and A being the symbol of finite differences with respect to n, it is plain that 

J"M n = ( 1) £k+n £ —(k+n+1)- 

Hence if the products £ k+n £' _(&+„+] , are computed for the various values 
of n, and are taken alternately with the positive and negative sign, and are 
written as if they were the successive differences of a function, we shall get 
the values of the factors M„ by filling out the scheme of differences. This 
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abbreviatiou is applicable equally whether we are making a numerical com- 
putation of the coefficient or a literal one. In the latter case the abbrevia- 
tion can be applied separately to each term of the form CeV" in the pro- 
ducts £"*£'_(,,+!). 

19. We proceed now to discuss the functions £. From their defini- 



tion we have 
whence 



k=z+ce / i \ 



(£)'cos> 



fc=+OS 

1 X' 



(-Jain /■ = *£_ 



£ 



f (i)- f (-D 



cos /<, 



sin /£. 



From which we gather that the functions £ can be computed by definite 

integrals thus £ ( j ) = — f [-)' cos(jv — kQdC. 
\ k J "• ° a 

Let us now suppose that the coefficient of «*, in the expansion of 
7]\¥(l — cosf- J ( I — — \ lJ 

in powers of s, is denoted by E I j J , then evidently 

\ J / fc=-» V i— ^ / ~2~ 

By writing in the expression 1 -t- s for s and changing the sign of j, it re- 
mains unaltered ; hence the relation 

By developing the factors of the expression 

7 V(1 _ a, s )W ( 1 — — J • 
by the binomial theorem, we easily get 

E ( i ) = (-1)Wp -i ' * -i] *' taM 

tx L 1+ r(^7+i) w + -xrn^T+iy^-i + 2) ; y +-J- 

This equation, as written, is correct only when k — j is not negative, but by 
h e relation given above we can reduce the case of Jc — j negative to that 
where it is positive. The factor in the brackets is a case of the series 
1+ aJ x+ x{a + 1) /?(/? + 1) g?+ 



l. r 



1.2. r .( r + l) 



— 1 77- 



treated by Gauss in a memoir entitled "Disquisitiones genei-ales circa seriem 
infinitum &o." (See Gauss' Werke, Vol. III. p. 123, and especially the 
"Nachlass." p.207). According to Gauss' notation 



<D- 



Whenever, of i -\- j and i — j, one is not negative, this series terminates af- 
ter a certain number of terms, thus affording a finite expression for the func- 
tion. But when these integers are both negative, the series is infinite. 
However it can be easily transformed into another which like the former is 
finite. From Gauss' investigation of these series, (See the volume just 
quoted, p. 209, Equation [82]), we have 

F(a, /?, r , x) = (\-x) a F( r - «, r — A r, 4 

Applying this to our expression, we get 



<t) 



E [i) = (-ir y p-y,A-j>/«^(i-^ 



2>+l 



XF(i + h + 1, * —j + l,k —j + 1, io 2 ). 
This expression is evidently finite when i—j and i -\-j are negative. 

20. The developments of the functions £ in powers of e. as far as ef 
have been tabulated by Prof. Cay ley in the Memoirs of the Royal Astro- 
nomical Society, A^ol. XXVII. It would conduce to the ready employ- 
ment of the preceding formulae if we had the function £ I j J explicitly 



expanded in ascending powers of e, but the attempts I have made to write 
such a series lead to extremely complex forms of the coefficients. Hence 
I shall give here only the coefficient of the lowest power of e in this func- 
tion, which suffices for obtaining all the terms of the lowest order in any 
coefficient of the expansion of 1 -+-J. We have when j — k is positive, 

and when h — j is positive, 

„.... +[i _, 0] ^](-r.' 
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21. Thus in the example alluded to above of the coefficient of cos (2^ 
4- 2/-), we find that the terms of the lowest order in f and £',(omit- 



\r> 



ting here, as in the scheme, the two lower indices), are 
£o = £ i = £2 = £3 = 1, 

[-2, 2]*+ [-2,1]^ 
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F_ r - [[-«, 3]-[-5, 2]f-f [-5, 1] ^--R, 0]^] (£ 
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48' 

590 
48' 

845 
48' 
1160 

IS"' 



Bringing into use our method of abbreviation, we multiply each of the 
preceding numerical coefficients by 48 in order to avoid fractions, and then 
write them alternately with the positive and negative signs in a diagonal 
line, and from these, as successive orders of 
standing in the vertical columns, thus, 
+ 389 



differences, derive the numbers 
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and dividing the numbers of the first column respectively by 1, — 1, 1 .2, 
— 1.2.3, we get the following as the terms of the lowest order in the coeffi- 
cient of cos (2£ — 5<T' -f 2f) in a'-i-J, 

1 I - (2) ,77,(2) 

i[ 3896 i + 2oifl i*+ 

which agrees with that found in the books. The following additional terms 
of the same coefficient can be written from the second, third, &c, columns, 
viz, those which are multiplied by e' 3 and the various powers of sin 2 J.7, 

^ -77H OVVa M ' -± 'JiOOa" -~4 -1- -' '«"■-"* ■'> -i- ff -'-:■ w>-mn" 

2 48 1 # da i 
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+ 2iii[ 845 ^ f +315a3 f t - 

— &c 
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5 ^(2,o)-i 
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-33a 1 ^„4 H- a 5 ^ § | e' 3 sin< 
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22. When we wish to obtain only the terms independent off and '(J , 
that is those on which the secular perturbations depend, i—0 and i'=0, and 
the Besselian function /disappears from the expressions giving the values of 

£ and £', and the coefficient of cos(i"y-\- i'"ff) in the expansion of -. can 

be written 

1 /i.'=ZO0 1IZZZ& f/Z-tl 

I I I ( 



i *?» -r- "?" , x ._„„+„_„ [&][«,j>] 



a' k—itnn=o j.=o ' w! 



X«' 



23. In leaving the subject of the development of 1-^z/, it may be well 
to note that two other forms can be given to the expression of the general 
coefficient, by employing, instead of the expression given above, either of the 
following 

e' I , , 1\ el , V 

--- — - # — a 



- /V-4 



a 



f("+7)-J('+i) 



,/(l_ (yV )(l- J 



But as they do not possess as much symmetry and brevity as the form giv- 
en above, we will pass over them. 

24. The second part of the Perturbative Function ,omitting the factor 
m' , is 

— ~ cos ip = — -U cos 2 ~. cos (v — v' + r) + sin 2 „ cos (v + v' + 6) 

— 2 r^ Sin 2L + A s J. 

According to the first theorem of § 3, the coefficient of z° in -s is 

equal to that of s° in 

^s(l — «w) (l— ^) 3 ; 

or it is equal to — 3rfco(l + w 2 ) = — |e. 

And according to the second theorem, the coefficient of z' in the same func- 
tion is equal to that of s' in 
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i4rwi-^) 2 ].^-i) 



■'^ &!-h 



i ds[_ 

=! (- 1Y 

Hence we have 

a t=-oo % L V y J 

And by simply writing 1-j-z for 2, 



r _, 0) /__l l=+00 



i=-o> *L 2 2—1 



The well known differential equations of elliptic motion 

supposing the axis of a; to be directed towards the perihelion, give us the 
equation 



«■*-. *(;•*) 



and consequently these two 
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^j^ — Ju y. 



By substituting these values in the expression given above for — —^camp 
it is not difficult to see that, in it, the coefficient of 

cos {ic + t'c + r) 

IS — COS 2 ^.^ J ie —uJie • »Y <" "'«' — ^iW 

and the coefficient of 

cos « + i'C + #) 



is r sin 2 - . 4 

a' 2 2 * 



C(i-l) (i+1) -| r- ii'-l) (i'+l)~| 

«/;« —C0V ie \.i'Tj' \ Jifet_ ~ (0 /2 J lJjL 



In the special case of i = the middle factors of these expressions take the 
indeterminate form 0-5-0, but then, in accordance with what has been show r n 
above, we should read — |e. Thus, by means of the Besselian functions, 
these coefficients take finite forms. 



